Project  Supervisor 
A. 


Department  of  the  Navy 
Office  of  Naval  Research 
Contract  N6onr-244.  Task  Order  H 


•*  1;  • '-y,  ^ .J 


'J.  ■ ‘ " 


* . \ • 


Supplement  to 


- - »v  -r  . ' * 0 . i . 


A.  J.  Acoeta 


- L . - 


’ *•'  ‘ v ;V 


• *•'  C‘<  *.  ’ y ' 


. • ‘ V - 

■ l 44  vV  >' 


<•  <?>:• . if, 


\ . A:>  <'.v. , -A 


4.  "■  ' ..  v . c''";  4-  l 4r’4444h 

.■*.  '•  - V"  ' ’ ' -X \ r.>  * ■ •’. •‘^V,v _ 4/  * 

~ t.  '*£&.$  • 4 -4)’  > 

. v?;..;-  • / '••i;V'v"  V^vTr  ,r.)  k ...,44  ■.>>“,  ..JT'  . < 6 *•  ’ „ 

>v" ■ * •>X  •*.  . / w * ' ■ ,v  *%>•>  ''*  • ''  *■*  *’;X'  J.  ^ - •-.»■•.  r:*. /• *n  .i  ' 

A ( .A  • Hy4«^yhani;c*  laborato.*  • 4>'J> 

c>;  V".  A.’'  .'.’j*®  • ''  *4.  >••*?•'*•  . "*  ''•  v-  • ’•■•!'  ' ' I*  * ••••  -f^r 

,Caljfwt;ht  Institute;  of  Technology  < X'i'  X ^ " A--;. 

■ ™ 5-  . , __v.  ..  _ ~,S . A A‘  ,r-»r  v. 

--  ' ''«c' ;^r ••  v.'r- ' ••'••'■  X;^15 


•'••■A”  •/.'•  '■••.'  % rA'-''"  • •'  ..'  ••'A'"1-;  ■ . 

' x''  . 4-^X.  fc-c,  • '*  ..'  * u > /I;'!  '.  . Vi  -e*  \ ’ '/./  ' -i'.K'. 

"S  > ' ~'-yt  —x  • *'  >**•  - •:*  ‘ "4-  ‘ . - S,  4 > ‘ H ^ if 

:.'  f V-;  V„.  ~ '/i:? 

A,  {■  Vr  ' a ;5-  V v v*  HolUtichfr  'vA!  -'W  ' 

! ' .•  : v < ••-•;■  ' . A : a ' • .'.  .-c-.-'-.v  -..A  • •/ 

,•■■  - :.  T v'V  l.vA  :•■•.  A'''.'  • 1 •“  'V' 

. ’ ■ -.  \.  . r*  "‘i  -•'•:•;•  •'••'•  " V;-./  ,-•••  : '•-A. 

:•  ' • >•  ' ‘ • *'/'>}  * X*K  ^ ‘.^y  r?rVn 

. Report.  ^ . • .V  ‘ •.  \>&y 

April.  1954  ' 'v  ';r.- 


*>  > , \ i 


v*  * . . 4' 


THIS  REPORT  HAS  BEEN  DELIMITED 
AND  CLEARED  FOR  PUBLIC  RELEASE 
UNDER  DOD  DIRECTIVE  5200,20  AND 
NO  RESTRICTIONS  ARE  IMPOSED  UPON 
ITS  USE  AND  DISCLOSURE, 

DISTRIBUTION  STATEMENT  A 

APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED, 


Abstract 


The  pffect  of  meridian  curvature  in  a centrifugal  impeller  is  estimated 
by  assuming  that  the  meridian  streamlines  may  be  replaced  by  "two- 
dimen  sional"  surfaces  of  revolution.  For  the  analysis,  a shape  was  chosen 
that  made  the  potential  flow  on  this  surface  easy  to  solve.  It  is  found  that 
the  slope  of  the  streamline  at  the  blade  inlet  is  quite  important  for  the  deter- 
mination of  the  flow-rate  for  shockless  entry.  If  the  solidity  is  somewhat 
greater  than  unity,  the  shockless  condition  is  relatively  independent  of  the 
exit  inclination.  The  head,  however,  shows  little  dependence  on  the  inlet, 
but  is  a function  of  the  exit  slope. 


Introduction 


The  theoretical  operating  characteristics  of  radial  flow  centrifugal  pump 
impellers  with  small  variations  in  breadth  and  vane  angle  have  been  worked 
out  previously.*  These  calculations  were  intended  primarily  for  straight 
radial  meridian  profiles  although  it  was  shown  how  they  could  be  applied  (in  the 
case  of  small  constant  breadth)  to  conical  shapes.  It  is  well  known  that  most 
impellers  of  the  so-called  Francis  type  have  appreciable  meridian  curvatures 
in  which  the  exit  is  more  or  less  radial  and  the  inlet  may  be  nearly  axial. 

This  present  note  is  a continuation  of  the  work  in  Reference  1 and  is  intended 
to  show  the  effect  of  this  curvature  on  the  characteristics  of  the  impeller. 

For  this  purpose,  the  meridian  flow  is  assumed  to  occur  on  a "two-dimensional" 
surface  of  revolution.  Furthermore,  two  adjacent  meridian  streamlines  are 
presumed  to  be  a constant  infinitesimal  distance  dh  apart  (Fig.  1),  so  that 
complex  variable  theory  can  be  used. 

If  dh  is  not  constant,  the  stream  function  v/ill  not  satisfy  L,aplace's 
equation  and  the  problem  becomes  much  more  difficult.  However,  in  this 
event,  further  corrections  can  still  be  made  as  are  indicated  in  (l).  The  im- 
peller may  be  thought  of  as  consisting  of  a number  of  such  surfaces,  so  that  a 
"psuedo  - three -dimensional"  solution  could  be  constructed  for  any  given  design. 
In  the  development  that  follows,  however,  only  one  meridian  stream  surface 
will  be  considered. 
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Formulation  of  the  Problem 

Consider  the  two-dimensional  surface  of  revolution  shown  in  Fig.  2 
where  the  equation  of  the  surface  is 


z - f(R)  . 


0) 


The  flow  on  thi3  surface  is  assumed  to  satisfy  all  of  the  requirements  of  a 

potential  flow,  i.e.,  it  is  frictionless,  ir rotational,  and  incompressible. 

Thus,  if  F is  the  potential  of  the  flow,  then  it  must  satisfy  Laplace's  equa- 

2 

tion  on  the  surface.  In  curv?.linear  coordinates  x^,  x^,,  this  equation  is 


L ( hA  *j[\  + _d_  ( W ; 

X1  V hl  Bxl/  dxz\h2  dx2  / 


0 . 


(2) 


Let  = 9 be  the  polar  angle,  and  x,  = s,  the  arc  distance  along  a meri- 
dian (0  = constant).  Since  an  element  of  distance  is 


2 2 2 2 
dS  = ds^  + R^de4 


then 


Thus 


h ■ = i ; h2  = R . 


R 


a2F 

57 


0 


where  R - R(s)  is  the  question  to  be  solved.  It  is  more  convenient  to  map 
the  (s,G)  coordinates  into  polar  coordinates  (y„  0),  where  Eq,  (2)  assumes 
the  form 


y 


JL 

ay 


o . 


O) 


The  connection  between  s,  R,  and  y is 


- J - 


or 


dy  d s 

y " R 


dR 

~R 


I + 


/ dz,  \ ^ 1 

idR'i  J 


1/2 


A1  30 


3 = 


(4) 


(5) 


The  relation  between  the  (s,0)  and  (y,  G)  planes  is  a conformal  one  so 
that  an  impeUer  blade  of  constant  angle  y defined  by 

tan  Y - R d9/ds 

is  mapped  into  a logarithmic  spiral  in  the  (y,  0)  plane. 


Boundary  Conditions 


In  the  (n,  0)  plane  let  V be  the  normal  velocity  component  and 
be  the  normal  component  in  the  (y,  9)  plane.  Then 


V 

ny 


or 


V 

ny 


V = V — — 

ns  oy  ns  3R  dy  ' 


= V 
ns 


R 

y 


(6a) 


V 

riy 


(6b) 
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For  an  impeller  the  total  flow  is  composed  of  the  through  flow  and  the 
displacement  flow.  The  through  flow  consists  of  the  flow  originating  from  a 
source  vortex  in  the  (y,  0)  plane  and  is  mapped  onto  the  physical  plane  by 
Eqs.  (4),  (5).  The  displacement  tiow  arises  from  rotation  of  the  vanes  with 
no  net  discharge.  At  any  point  on  the  blade  surface  the  normal  velocity  com- 
ponent must  satisfy  the  condition  (Fig.  2) 

V = Rocos  r . 
n s 


Thus, 


V 

ny 


(7) 


is  the  boundary  condition  imposed  in  the  (y,  0)  plane.  In  this  plane,  the 
blade  will  appear  as  a logarithmic  spiral  for  which  suitable  transformations 
are  available. 


Applications 


A.  Francis-type  Impeller 

Ir  this  type  of  impeller  the  exit  is  nearly  radial  and  the  inlet  may  vary 
anywl  ie  from  axial  to  near  radial  as  seen  in  Fig.  1.  From  Eqs.  (4)  and 
(7)  it  is  seen  that  rather  difficult  boundary  conditions  may  occur  in  the  (y,  0) 
piano  depending  on  the  exact  shape  of  the  streamline  contour.  In  this  work 
the  item  of  major  interest  is  the  axial  inclination  of  the  inlet  (®j)  and  not  the 
particular  shape  by  which  this  slope  is  achieved.  It  is  supposed  for  the 
present  that  the  details  of  the  profile  shape  are  not  of  crucial  importance, 
particularly  for  impellers  of  high  solidity.  Later,  this  contention  will  be 
supported  by  comparing  the  results  of  two  different  shapes. 

In  view  of  these  remarks  it  seems  reasonable  to  choose  that  contour 
which  will  give  the  simple  st  boundary  condition  in  the  (y,  0)  plane  and  still 
give  a streamline  shape  that  has  a continuous  smooth  curvature  from  inlet  to 
exit.  A simple  relation  that  still  preserves  the  axial  to  radial  transition  is 
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c. 

y 


where  R.  is  a constant.  Equation  (7)  gives 


V = to  co  s Y 

ny 


(8) 


(9) 


The  first  term  is  the  boundary  condition  for  an  axial  impeller  (or  cascade) 
of  straight  blades.  The  second  term  is  the  boundary  condition  for  a strictly 
radial  impeller. 

The  relation  between  z and  R is  givenby  Eq.  (4),  i.  e.  , 


dX  = 
y 


dR 

R 


. 


1 + f 


*1 


1/2 


or  with  Eq.  (8) 


f’(R) 


dz 

dR  = 


R 


L (RZ-Rf)2 


i/2 


(10) 


The  integral  of  Eq.  (10)  is  plotted  in  Fig.  3 and  the  angle  is  given  in 

Fig.  4 as  a function  of  R/R.  . The  constant,  R.,  is  the  value  of  the  radius 
for  which  the  streamline  becomes  axial.  Although  the  shape  given  in  Fig.  3 
has  a much  smaller  radius  ratio  than  is  typical  of  pump  designs,  it  will  be 
seen  that  it  is  still  useful  for  the  determination  of  the  effect  of  6^  . 


1.  Displacement  Flow 

Inorder  to  present  the  results  in  dimensionless  form,  consider  the  r.cw 
variable 


Then 


V 

nr 


V 

= 1 

Rz 


d(y/R2) 


dr 


(12) 
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or 


With  these  substitutions  r = 1 when-  R/R_  = I,  and  all  velocities  are  now 

£. 

expressed  in  terms  of  the  tip  speed  R of  the  impeller. 

SLm 

The  blades  appear  as  logarithmic  spirals  in  the  (r,  0)  plane,  and 
because  of  that  fact  this  plane  can  be  easily  transformed  into  one  in  which  the 
blade  system  appears  as  a circle.  In  this  latter  plane  the  solution  to  the 
boundary  value  problems  given  by  Eq.  (13)  is  relatively  simple.  As  pre- 
viously noted,  the  second  term  represents  the  conventional  radial  flow  theory 
whereas  the  first  term  is  the  same  as  that  due  to  a source  flow  at  the  origin. 
Both  of  these  results  are  given  ir.  Kef.  I in  a form  most  useful  for  presenta- 
tion here. 

In  this  way  the  effect  of  meridian  curvature  can  be  obtained  by  ''fairing 
together'*  two  known  solutions  in  the  (r,  8}  plane  by  means  of  Eq.  (13).  The 
relation  between  the  actual  geometry  and  the  transformed  radial  plane  (r,  0) 
is  given  by  Eqs.  (4),  (5),  and  (11) 

The  solution  of  the  second  term  is  given  by  Eq.  ( 1 I ) in  Ref.  2.  If  sub- 
script {2)  denotes  the  exit  and  (i)  the  inlet,  then  the  tangential  velocity  in  the 
circle  plane  due  to  this  term  is 


where  w , p are  constants  that  occur  in  the  transformation  of  the  (r„  ft"5  into 
the  circle  plane.  The  constant  q is  a function  of  the  number  of  blades  N 


and  the  blade  angle  y and  F.,  F,  are  hyper  geometric  functions  given  by 
Eqs.  (13)  in  Ref.  !. 

The  tangential  velocity  of  the  first  (or  source  flow)  term  is  readily  ob- 
tained in  the  circle  plane  by  letting  m = -2  or  N*  = oo  an  Eq.  {ibj.  Ref.  1, 
to  obtain 


2 cosy 
* P, 


_ 2 cosY  i 

'*Z  N Pi  L r2- 


(15) 


The  dashes  refer  to  the  first  term  of  the  displacement  flow  given  in  Eq.  (13). 


2.  Through  Flow 

The  through  flow  solution  in  the  (r,  ©)  plane  is  not  altered.  In  parti- 
cular, the  flow-rate  coefficient  4 - V /V^  remains  unchanged  in  going  from 

r 2,  “ 

the  physical  plane  to  the  (r,  ©)  plane.  The  appropriate  solutions  are  then 
given  by  Eqs.  (15),  Ref.  I,  i.  e.. 


. Q ( tan  Y - tan  a ) 

n p2  1 

(16) 

T~-J'anT  - tanaj) 

where  Q is  the  flow- rate  per  passage  per  unit  width. 

3.  Shockless  Entry 

All  the  essential  information  is  now  available  to  make  calculations  of 
head  flow-rate  performance  in  which  the  effect  of  meridian  curvature  is  taken 
into  account.  The  first  and  most  important  item  is  the  flow- rate  for  which 
stockless  or  smooth  entry  occurs  on  the  Francis-type  impellers. 

This  operating  point  is  determined  by  the  condition  that  the  velocity  must 
be  zero  in  the  circle  plane  at  the  point  corresponding  to  the  inlet  of  the  blade. 
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According  to  Ref.  I this  sc^ui  s when 


(V.  - V'  ) * {V.  - V.  ) + (V,  -V  i = 0 

‘l  >2  dl  d2  dl  “2 


(17) 


After  suitable  manipulation  of  the  equations  quoted  above,  the  flow- rate  co- 
efficient 4 for  this  condition  may  be  expressed  as 


where  $ is  defined  by 


V 

r2  _ Flow-Rate 
U2  " Area  x Tip  Speed 


(18/ 


(19) 


TSie  subscript  (oo)  denotes  the  value  of  this  coefficient  for  an  infinite  number 
of  vanes,  i.  e. , 


4 


e ao 


- tan  Uj)  - 


The  term  A is  the  ratio  4 /<$  for  a purely  radial  impeller,  and  is  given 

e e oo 

in  Ref.  1 as 


A 


0 


T 


) 


(20) 


subject  to  the  restriction  that  the  solidity  be  greater  tlt&u  about  1.2.  The 
solidity  for  a Francis  impeller  is  defined  as 


N In  r Tj 
2x  cos  y 


(21) 


Equation  (18)  appears  to  be  independent  of  Rj/R^  . Actually,  A depends  on 


i lard  thus  H P hut  if  I he  solidify  is  subjected  to  the  above  restriction 
-4  becomes  bndepe?iflcnl  of  'jr  and  is  a function  of  the.  vane  angle  and  number  of 
blades  alone.  In  fact,  Eq.  (20)  is  only  an  aTjproxirratioin  for  the  case  sr  = i.  2, 
and  far  smaller  values  a more  complicated  relation  must  be  used  (Ref.  1). 

If  a large  4'alue  of  F_  R.  is  chosen  as  the  tip  radius  in  Fig.  3,  the 
angular  inclination  of  the  exit  will  be  essentially  radial.  Suppose  for  the  mo- 
ment that  the  radius  raho  R R7  of  the  impeller  is  sufficiently  small  so  that 

t c 

*r  — 1.2.  If  now,  R R,  is  decreased  below  this  value,  the  inlet  edge  of  the 
blade  will  become  more  axial.  Since  the  impeller  geometry  already  i--  pre- 
sumed tc  be  such  that  A does  not  change,  the  sole  effect  on  the  shockless  flow- 
rate  condition  is  that  due  to  a change  in  the  inlet  meridian  angle  6 . Values  c-f 
(R.  -'P.)  then  may  be  substituted  into  Eq.  (13)  and  the  corresponding  value  of 
5 ^ determined  from  Fig.  4.  (Actually,  for  these  purposes,  one  may  just  as 
well  take  R_/R.  as  infinity.  ) 

This  result  is  independent  of  the  exit,  meridian  angle  G if  the  solidity 
remains  high.  This  fact  may  be  seen  in  the  following  way:  Let  <r  have  a 

constant  value  so  that  Tj  - r j/r7  is  also  constant.  Since 


the  value  of  A will  stay  constant  if  Tj  is  kept  constant  and  independent  of  the 

radius  ratio  .R  j R^  of  the  impeller.  Thus,  for  a given  r j (and  therefore 

constant  A)  R../R  and  R /R  can  both  approach  unity  and  G_  approach  90° 
l IS  lb  I L 

without  affecting  the  value  of  ^ as  computed  by  Eq.  (18). 

Equation  (18)  is  -ilotted  in  Fig.  5 for  the  case  of  N - 6,  Y =70°  which  is 
typical  of  many  pump  designs.  It  i s evident  that  the  inlet  meridian  inclination 
fi  has  a large  effect  and  must  be  taken  in  account  in  any  actual  design. 

4..  Head  Flow-Race  Performance 

The  head  developed  by  the  impeher  is  usually  expressed  by  the  coef- 


ficient 
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where  H is  the  developed  head,  and  is  the  tip  speed  of  the  impeller. 

This  coefficient  may  he  shown  to  be  equal  to 


♦ 


co  and  r _ being  unity.  With  the  aid  of  Eqs.  (14),  (15),  and  ( 16),  the  head 
flow-rate  relation  becomes 


| - * CHfi(tanY  - tan  c ) (22) 

for  backward  curved  impellers  where 


is  now  the  shui-off  head  coefficient.  The  coefficient  depends  on  the  im- 

peller solidity  and  the  blade  angle  and  is  usually  very  near  unity.  The  term 
f is  the  shut-off  head  coefficient  for  a straight  radial  configuration.  Graphs 
of  both  of  these  coefficients  are  given  in  Ref.  1. 

Equation  (23)  is  plotted  in  Fig.  7 for  N = 6 and  y = 70°  with  the  assump- 
tion that  et  = 1.  2. 

B.  Conical  Impellers 

The  meridian  picture  now  appears  as  a simple  cone.  Equations  (l)  and 
(4)  give 


F(R)  = tan  0 


- \ J_]COST 
’ i RZJ 


The  boundary  condition  for  the  displacement  flow  is 


V_ ^ = » cosy  r 


2 cos  5-1 


(24) 


The  solution  of  this  problem  is  given  by  a value  of 


m = Z cos  6 - 2 

in  Eq.  (35),  Ref.  1.  It  is  seen  that  the  equivalent  number  of  vanes  to  be  used 
in  Eq  ^38),  Ref.  1,  is 


N'  = 


2 N 


(25) 


Thus,  the  characteristics  of  the  conical  impeller  are  obtained  by  replacing  the 
actual  number  of  blades  by  the  equivalent  number  given  by  Eq.  (25).  The 
shockless  entry  calculation  gives  then 


e oo 


TJ n7! 


4tt 


o s 

F 


3 N 


(26) 


These  results  were  previously  obtained  in  Ref.  1 by  a different  method  and 
also  by  Prof.  Wislicenus  in  an  unpublished  note. 

Equation  (26)  is  plotted  on  Fig.  5 for  the  case  of  N = 6,  y * 70°.  It 
is  seen  that  this  curve  is  quite  close  to  that  of  Eq.  (18).  This  agreement 
shows  that  as  far  as  shocklcss  entry  is  concerned  the  exact  details  of  the 
meridian  streamline  shape  are  relatively  unimportant  and  that  the  influence 
of  the  exit  angle  is  small. 

The  shut-off  head  coefficient  for  a conical  impeller  is  determined  by 
the  use  of  the  equivalent  number  of  blades  as  given  by  Eq.  24  and  Fig.  6 of 
Ref.  1.  The  effect  of  the  cone  angle  on  this  coefficient  is  shown  in  Fig.  6 
for  the  case  N = 6,  Y = 70°  and  <r  = 1.2.  The  two  curves  are  seen  to  agree 
quite  closely  as  is  reasonable  tor  impellers  of  nigh  solidity. 

Conclusions 

This  analysis  has  shown  that  the  fiow-rate  for  shockless  entry  is 
strongly  dependent  upon  the  inlet  meridional  inclination  of  the  blade  and  is 
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relatively  independent  of  the  shape  of  the  me ridionat  stream  surface.  The 
ahut-off  head  coefficient  shows  a similar  behavior  with  the  exit  slope. 

These  calculations  were  made  for  a particular  profile  that  hid  ? con- 
tinuous axial  to  radial  curvature,  and  for  conical  profiles.  The  agreement 
between  the  two  was  sufficiently  close  so  that  the  effect  of  inlet  and  exit 
meridional  slopes  may  just  as  well  be  approximated  by  the  simpler  conical 
shape.  This  result  should  hold  true  at  least  for  shroud  curvatures  no  sharper 
than  the  particular  one  studied, 
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Notations 

A - The  flow-rate  ratio  </>  / d for  a radial  impeller 

e e oo 

- Flow-rate  correction  factor 
F - Velocity  potential 

g - Gravitational  constant 

H - Impeller  head  (ft.  ) 

N - Number  of  blades 

•3 

Q - Flow  rate  (ft  /sec) 

R Radial  coordinate  in  physical  plane 

r - Radius  in  the  transformed  radial  plane 
s - Arc  distance  along  a meridian 

U - Circumferential  velocity'  (R&>) 

V - Absolute  velocity 

y - Intermediate  radial  variable 

z - Axial  coordinate  in  physical  plane 


a 

Y 

3 
9 

a 

4 


0) 


- Angle  between  absolute  velocity  and  meridian  (positive  increasing 
counter-clockwi  se) 


- Angle  between  blade  tangent  and  meridian 

- Angle  between  radial  line  and  tangent  to  meridian 


- Polar  angle  about  axis  of  rotation 
N In  r^ 


- Solidity  - -= — 
’ Z IT 


cos  y 

Flow-rate  coefficient. 


Di  schar ge 


Exit  area  x tip  speed 


- Head  coefficient,  H/U-^/g 

- Shut-off  head  coefficient 

- Shut-off  head  coefficient  for  a purely  radial  impeller 

- Angular  speed 


Refers  to  displacement  flow 
Denotes  shockless  or  smooth  entry 
Normal  component  of  velocity 
Meridional  component  of  velocity 
Refers  to  (r,  0)  plane 
Refers  to  (s,  0)  plane 
Refers  to  through  flow 
Refers  to  vane  inlet 
Refers  to  vane  exit  or  tip 
Denotes  the  infinite  vane  theory 


Fig.  1 - Meridian  profile  of  a Francis  type  impeller. 


Fig.  2 - Oblique  view  of  a typical  stream  surface  of  revolution 
showing  the  blade  trace  and  coordinate  system. 


ANGLE  BETWEEN  MERIDIAN  AND  RADIUS 


Fig.  3 - Plot  of  axial  distance  vs.  radius  for  the  special 
stream  surface  given  by  Eq.  10. 
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Fig.  4 - Meridian  inclination  <5  for  the  shape  shown  in  Fig.  3. 


SHUT-OFF  HEAD  COEFFICIENT,  • RATIO  OF  SHOCKLESS  FLOW  RATE 
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5 - Ratio  of  shockless  flow  rates  vs.  inlet  meridian  shroud  angle  6j 
for  the  case  of  6 blades,  Y = 70°  anc*  solidity  =1.2. 


Fig.  6 - Shut-off  head  coefficient  <jr  vs.  exit  meridian  shroud  angle  5, 
for  the  case  of  6 blades,  Y *=  70°  and  solidity  = 1. 2, 
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